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A cross-power spectrum phase based adaptive technique is discussed which iteratively 
determines the time delay between two digitized signals that are coherent. The adaptive 
delay algorithm belongs to a class of algorithms that identifies a minimum of a pattern 
matching function. The algorithm uses a gradient technique to find the value of the adaptive 
delay that minimizes a cost function based in part on the slope of a linear function that 
fits the measured cross power spectrum phase and in part on the standard error of the 
curve fit. This procedure is applied to data from a Honeywell TECH977 static-engine 
test. Data was obtained using a combustor probe, two turbine exit probes, and far-field 
microphones. Signals from this instrumentation are used estimate the post- combustion 
residence time in the combustor. Comparison with previous studies of the post- combustion 
residence time validates this approach. In addition, the procedure removes the bias due 
to misalignment of signals in the calculation of coherence which is a first step in applying 
array processing methods to the magnitude squared coherence data. The procedure also 
provides an estimate of the cross-spectrum phase-offset. 


Nomenclature 

a Zero intercept of two parameter phase function 4>' k m (f) 

B e Bandwidth resolution, B e — A / = l/Td = r s /N, 16 Hz 
b k m Slope of cross-spectrum phase curve 

C Q , CYiax, C'min, C Cost functions 

Dp m Propagation time from turbine exit to far-field microphone 
D k t Post combustion residence time in the combustor, D k t = D k m — -Df m 
D k m Propagation from combustor to far-field microphone 
E [ ] Effective Value of E [ ] 

/ Frequency, Hz 

f u Upper frequency limit, f u = l/2At = r s / 2, Hz (32 768 Hz) 

F 2 Weight function 

Gee(f) Estimated output signal one-sided auto spectral density function 
G k k(f) Estimated input signal one-sided auto spectral density function 
Ckm(f), Ctm(f) Estimated one-sided cross-spectral density function 
J Cost function 

N Segment length, number of samples per segment, N = r g T^. (4 096) 
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n Signal noise 

nk£ Number of points used in two parameter curve fit 

n D Number of segments/blocks with 50 percent overlap, n 0 — 2 S e T tota i ~ 2240 

n s Number of disjoint (independent) segments used in spectra estimates, n s — B c T tota \ ~ 1120 

r(t ) Acoustic signal 

•Rfcm(T') Cross-correlation function 

r s Sample rate, samples/s (65536) 

S Cost function sums 

s(t) Acoustic signal 

S j ] Summation function over [ 

t Time variable 

Ttotai Total record length, T tota i = n s Td, s (ss 70 s) 

Td Record length of segment, Td — N/r s , 0.0625 s 
Wk Weight factor 

Xj Frequency fj 

yj Measured cross-spectrum phase 0™ m 

Subscripts 

1 Combustion location 

2 Turbine exit location 

3 microphone location 

£ Location is turbine 

k Location is combustor 

m Location is far-field microphone 

x Input signal 

y Output signal 

Symbols 

Mean value of [ ]. 

< > Angular brackets denote long time averages. 

Sk m Slope of two parameter phase function <f/ k m (f) 

At Sampling interval, At = 1 fr a s 

e Residual 

ll.y(f) Estimated magnitude-squared coherence (MSC) function 
Ta„,i/ n ( n s) MSC threshold function 
Estimate of [ ] 

S Imaginary part 

v weight value 

u) Angular frequency, u> — 2 n f 

4> Angle variable 

4>' k m (f) Two parameter linear curve fit function 

0o Zero intercept of two parameter phase function (j)' k m (f) 

0fc m{f) Cross-spectrum phase function 

%b Weight factor 

9? Real part 

<7 k m Residual variance 

at m Standard error 

r Time displacement 

Int Integer part 

m Measured 


I. Introduction 

The technique introduced in this paper is an adaptive estimation process which synchronizes two 
signals correlated over a frequency range using a time domain iterative algorithm. This procedure will enable 
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removal of the bias due to misalignment of signals in the calculation of coherence and will estimate more 
accurately the cross-power spectrum phase-offset. 

The algorithm uses a cost function and a gradient approximation to continually update the time 
delay estimate until the process has minimized the cost function. The cost function is based on a linear 
curve fit to the unwrapped phase of the cross-spectrum between the two signals over the frequency range of 
interest. The linear curve fit assumes the existence of a cross-spectrum phase-offset and does not assume the 
linear fit has a zero intercept at zero frequency. The iterative technique used is a steepest-descent gradient 
technique and thus requires that the cost function is unimodal. 

A. Relevance 

Low frequency noise generated in the turbofan engine core can make a significant contribution to the overall 
noise signature in the aft direction at the low power settings which are used on an airport flight approach 
trajectory. This type of low frequency noise may become even more of a problem for future aircraft. Two 
possible low frequency noise sources are direct and indirect combustion noise. The source of combustion noise 
attributed to the unsteady pressures produced by the unsteady combustion process that propagate through 
the turbine to the far field is called the direct combustion noise source. The other source of turbofan engine 
combustion noise is known as the indirect mechanism in which the noise is generated in the turbine by 
the interaction of entropy fluctuations, which also originate from the unsteady combustion process, as they 
propagate through regions characterized by mean flow velocity or pressure gradients in the turbine stages. 
This indirect noise source was studied using analvtical models by Ffowcs Williams and Howe, 1 Pickett, 2 
Marble and Candel, 3 * Cumpstv and Marble, 4,5 Cumpsty, 6 Gliebe et ah, 7 Mani, 8 Bodony 9 and Leyko 10 . 

The net travel time of the indirect combustion noise signal from the combustor to the turbine exit 
and the far field is longer than the direct one since the travel velocity of the entropv fluctuations to the 
turbine is the flow velocity in the combustor. This flow velocity is a small fraction of the speed of sound. 
Miles et al. 11,12 has shown the pressure and entropy should be in phase in the combustor. Consequently, one 
might expect that the pressure signal from an indirect combustion noise source would be delayed relative to 
a pressure signal from a direct combustion noise source since an indirect combustion noise signal does not 
travel with the speed of an acoustic wave until it interacts with the turbine. Miles 13-15 with data from the 
Honeywell TECH977 engine test program 16 shows that the cross-spectra and correlation function between 
a combustor sensor and far-field microphones are tools that provides a way to separate low frequency direct 
and indirect coherent combustion noise due to this travel delay time. 

B. Adaptive Estimation 

Etter and Stearns 17 describe an autocorrelation function based adaptive estimation process for time delays 
for possible use with sonar or radar signals. A similar procedure is outlined by Carter 18 for passive sonar 
signal processing which uses a cross-correlation function. The generalized correlation method for estimation 
of time delay that uses a pair of prefilters before a cross-correlation process has been discussed by Knapp 
and Carter 19 and Azaria and Hertz. 20 

These methods are early versions of modern signal-processing algorithms now widely used in medical 
ultrasound, speech processing and other applications that can be more generally described as the identifi- 
cation of a minimum (or maximum, depending on the particular algorithm used) of a pattern matching 
function. 21-23 A great many pattern matching algorithms are used. Many are reviewed by Jacovitti and 
Scarano 21 and Viola and Walker. 23 While some cost functions use cross correlation, in many applications 
other signal matching measures have been developed for applications requiring intensive or real-time pro- 
cessing. 

C. Present Technique 

The adaptive estimation process used herein estimates the time delay with a cost function that is based 
on the unwrapped cross-spectrum phase between two sensors over a preselected frequency range. The cost 
function is based on the slope of the regression line and the standard deviation of the linear curve fit to the 
phase. It is assumed that the time series overlap for a time delay that minimizes this cost function. The 
procedure at each iterative step determines the value of the cost function and continually realigns the time 
series by varying the time lag until the cost function has been minimized and thereby effectively has reduced 
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the slope to zero. The frequency inter\al is preselected so that the magnitude-squared coherence between 
signals is not too small in the selected range. A non-adaptive time delay estimation method that uses a 
linear fit with a zero intercept at zero frequency to unwrap cross-spectrum phase is described by Piersol. 24 
In contrast, the method used herein does not require the measured phase shift sequence to pass through the 
origin which is an assumption used by Piersol 24 and other investigators. 25 27 Since the method uses only 
a selected portion of the calculated phase, it is called a generalized unwrapped phase method. The same 
results could be achieved at more computational expense by using a pair of pre-filters to select the frequency 
range of interest as is done by Knapp and Carter. 19 Other non-adaptive methods that estimate a time delay 
using the cross-spectrum phase are available. 25-27 

This adaptive estimation method of calculating the time delay is applied to time series measurements 
in the combustor, at the turbine exit, and at select far-field microphones from the Honevwell TECH977 test 
program. 16 In order to calculate the post-combustion residence time in the combustor i.e. the time delay 
of the indirect combustion noise in the combustor, the propagation time from the combustor probe and the 
turbine exit probe to a far-field microphone are directly estimated by the procedure and subtracted. This 
is the first application of time delay estimation using a phase based adaptive pattern matching function 
to measure the post-combustion (post-flame) residence time in a turbofan engine combustor using signals 
measured at the turbine exit with signals from a combustor sensor probe and far field microphones. 

II. Engine test data 

The NASA/Honeywell static engine test program was conducted at Honeywell’s San Tan outdoor 
acoustic test facility using a Honevwell TECH977 engine (Fig. 1) and the results are described in a report 
by Weir. 16 The dual-spool, turbofan engine has a direct drive, wide chord fan connected by a long shaft to 
the low-pressure turbine spool and a high-pressure compressor connected by a concentric short shaft to the 
turbine high-pressure spool. The fan diameter is about 0.87 m. The combustor design is a straight-through- 
flow annular geometry with 16 fuel nozzles and 2 igniters. Data obtained for one configuration in the test 
program is analyzed in this paper. The engine-internal instrumentation in this configuration included high- 
temperature pressure sensors with air cooling in a combustor igniter port identified herein as CIP1 and at 
the turbine exit sensors identified as T551 and T552. Pressure time histories from these internal sensors and 
far-field microphones are used herein. 

The data acquisition system had a sampling rate of 65 536 Hz and a sampling duration of roughly 
70 s. The spectra were calculated using a 50 percent overlap. This permitted data reduction using ~ 1120 
non-overlapped ensemble averages or with a 50 percent overlap « 2240 overlapped ensemble averages at a 
bandwidth resolution of 16 Hz. Further signal estimation parameters are shown in Table 1. 

The engine condition power settings discussed are 48 and 54 percent of the Max Power setting. 
Microphone locations used are at 90°, 110°, 130°, and 160° measured from the inlet. The test was conducted 
when the air temperature was about 9° C (48.2° F) and the microphone radius was 30.48 m. The engine 
was at a height of 3.049 m. 

As discussed, 13,28 the size of the magnitude-squared coherence between the combustor sensor and 
a far-field microphone identifies the spectral region of importance for indirect combustion noise as being in 
the 0-200 Hz frequency range. Consequently, herein the cross-spectral phase measurements between engine- 
internal sensors and a far-field microphone only over the frequency range 16-208 Hz is used to estimate the 
propagation time to a far field microphone. 


III. Theory 


A. Procedure 

The far-field microphones receives noise from the turbofan, jet, and core. The procedure characterizes the 
core noise as indirect combustor noise and turbine-exit coherent broadband noise. Furthermore, the method 
estimates the propagation time from source to far-field microphone. A simple mathematical model for the 
propagation time measurements is 

rk{t) = si(t — Dk m) + ni(t) Combustor Sensor (1) 

ri(t) = S 2 {t — Dg m ) + ri 2 {t) Turbine Sensor (2) 
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r m (t) — s 3 (t) + n 3 {t) Microphone 


(3) 


where the signal at station k from the combustion pressure sensor is delayed by the propagation time delay, 
D km , the turbine pressure sensor signal at location i is delayed by the propagation time delay, Di m , the 
signal at station m is from the far-field microphone. 

The signal at station k is from the combustor and the acoustic pressure signal s\{t) is assumed to 
be related to an entropy or temperature disturbance. This disturbance moves at the combustor flow velocity 
to the turbine where it interacts with the pressure field creating indirect combustion noise which reaches 
a microphone after a total travel time, D km . The signal at station C is from the turbine. The acoustic 
pressure signal S2(t) is from the turbine exit coherent broadband noise which reaches a microphone after a 
total travel time, D f rn . In either case the procedure to finding the time delay is identical and is described in 
the algorithm section. The signals at station k, £ and m are contaminated by independent, Gaussian, and 
stationary noise terms, rai(f), 712 (t), and 773 (f), which are uncorrelated with si(t), S 2 (f), and 53 (f) and each 
other. 

The cross-correlation function between the combustor signal and a far-field microphone is given by 

Rkm(r) = E[r k (t)r m (t + r)] = R SlS3 {r - D k m ) (4) 

where R SiS3 (t — D k m ) is the autocorrelation function with a peak at t = D k m . 

In the frequency domain, the Fourier transform of equation (4) is 

Gkm(f) = G SlS3 (/)exp(- i 2tt / D km ) (5) 

where G SlS3 is the one-sided cross-spectrum and the time delay D km appears in the cross spectrum as a 
phase function 


$k m(/) — 2 7T / D k m — U) D k m 

where oj = 2 n f . The magnitude-squared coherence (MSC) is defined by 

|Gfcm(/)| 2 


7 km(f) = 


(6) 


(7) 


by 


|G fcJfc (/)||G mm (/)| 

The cross-correlation function between the turbine probe signal and a far-field microphone is given 


R(. m{r) = E[re(t)r m (t + r)] = H, 2S3 (r - D t m ) (8) 

where R S2S3 (t — Dt m ) is the autocorrelation function with a peak at t — Dt m . 

In the frequency domain, the Fourier transform of equation (8) is 

Gtm(f) = G S2S3 (f) exp(— i2n f Dtm) (9) 

where G Sl S3 is the one-sided cross-spectrum and the time delay Dt m appears in the cross spectrum as a 
phase function. 


4>tm(f) =2 7T / D tm =u Dtm (10) 

The magnitude-squared coherence (MSC) is defined by 


7 tm(f) 


\Gtm(f )\ 2 
\Gee{f)\ \G mm (f)\ 


( 11 ) 


Then D k t — D k m — Dt m . With a small error since the acoustic propagation speed are high and the distances 
involved are small, the indirect combustion noise is dealt with as if it all begins propagating acoustically at 
the turbine exit and then it propagates to a microphone in the same manner as the direct combustion noise. 
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B. Signal Processing Method 

All the spectra and cross-spectra are estimated using Welch’s nonparametric method which is based on 
averaging multiple windowed periodograms using overlapping time sequences. 29-33 Using these spectra 
and cross spectra, the magnitude-squared coherence (MSC) is calculated to measure the similarity of the 
amplitude variations at particular frequencies. 33 In some sections, the “ accent will be used to denote the 
statistical basis of a variable and that it is based on a finite time series. 


C. Accuracy of Estimates 

We consider the cross spectrum magnitude estimate, Gk m(f), and phase estimate, <frk m(f), computed from 
Fourier transforms of signals ri(f) and r^t) calculated for a sample length T tot ai divided into n s independent 
samples using a periodogram procedure. The estimated cross-spectrum Gkm is a complex number. The 
estimated cross-spectrum phase (f>k m is given by 


4>k m 


— arctan 


m 

$lG k m 


( 12 ) 


The spectral estimates are statistically independent of one another and have a standard deviation 
discussed by Piersol 24 , derived by Bendat, 34 and also presented in a book by Bendat and Piersol 35 approx- 
imately given by 






j \Gkm{fi)\ 
V n S ikmifi) 


ircsm 




(fA 


1/2 


(13) 

(14) 


For small phase errors when sin o = a 


&[&k m(fi )] 


' 1 ~ 7 ImUi) ' 
. 2n s 7 kmifi). 


1/2 


(15) 


It can be seen from Eqs (14)-(15) that the error in the estimated phase angle becomes infinitely large 
as the MSC tends to zero. This is the reason why only a limited frequency interval is used in the generalized 
adaptive phase method. 


D. True MSC based on the Adaptive Estimation Procedure 


To a first order approximation the computed/estimated unaligned and aligned MSC based on the far- field 
microphones and probes in the combustor and at the turbine exit are related by 


-2 
Ik m 


I unaligned 


-2 
It m 


I un aligned 


1 - 


D k 


j< d J 7fc m (/ ) I aligned 


1 - 


D t 


J h m(f) (aligned 


(16) 

(17) 


if Dk rn and Dz m are less than (Bendat 34 and Bedat and Piersol 35 ). 

To increase the number of blocks of data sampled to get better statistics the bandwidth resolution, 
B e , used is 16 Hz and the segment length, of the signal is then 62.5 ms. However, the time delay from 
the engine probes to the far field-microphone, Dk m and D? m is between 70 m s and 90 ms. Consequently, 
the alignment procedure must be used since using non-aligned measurements one would only measure the 
MSC threshold which is small. An order of magnitude estimate of the threshold for a 95 percent confidence 
interval and n s = 1120 independent blocks of data is (Miles 36,37 ) 


7aL,2/n ( n s) = 1 - (1 - Q.95) 1 /^- 1 ) = 2.67 X 10 -3 . 


( 18 ) 
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By aligning the measurements and removing any time delay the true MSC can be measured. The 
generalized unwrapped phase adaptive estimation procedure does need an initial condition for D k m and D? m . 
An estimated value of 70 m s to 90 ms based on the acoustic speed of sound and the 100 ft microphone 
radius is sufficient for this purpose. 


E. Algorithm 

The current algorithm identifies the dominant noise source at the turbine exit as either direct or indirect 
combustion noise over a certain frequency range. It does this by calculating the difference in delay time, 
Dkii for a signal from the combustor to reach a far-field microphone and a signal from the turbine exit to 
reach the same far-held microphone. If both time delays are about the same, the dominant noise source is 
direct combustion noise. If the time delay difference is in the region of 3 to 4 ms the dominant noise source 
is indirect combustion noise. 

The general idea is that the noise from the turbine is acoustic while the noise in the combustor has 
a spectrum related to the noise created by entropy waves interacting with the pressure held in the turbine. 
This interaction produces noise after a time delay, D ki , related to the convection velocity in the combustor 
and the distance from the flame zone to the turbine. This post-combustion (post-flame) time delay, D k e, 
shows up in the phase of the cross-spectra, of the combustion noise signal and a far-held microphone signal. 

The procedure used is to hrst evaluate the time delay from the combustor and far-held microphone 
cross spectrum, D k m . This is done by iteratively and optimally resetting the origin of the time history of 
the far-held microphone signal and at each step calculating the pressure cross-spectrum phase between the 
combustor probe signal and the far held microphone signal until the slope of the phase plot curve is zero. This 
time delav is the acoustic travel time, D^ m , plus the post-combustion (post-flame) time delay, D k g. Next 
the same method is used to evaluate the time delay from the turbine exit probe signal and the same far-held 
microphone, Dg m . Again this is done by iteratively and optimally resetting the origin of the time history of 
the far-held microphone signal and at each step calculating the pressure cross-spectrum phase between the 
turbine exit probe signal and the far held microphone signal until the slope of the phase plot curve is zero. 
This provides an estimate of the acoustic travel time from the turbine exit to the same microphone, Dg rn 
Consequently, the difference in these actual time delays is the post-combustion (post-hame) time delay, D ke . 
The indirect combustion noise is treated as if it all begins propagating acoustically at the turbine exit and 
then it propagates to each microphone in the same manner as the direct combustion noise. The resolution 
is limited bv the sampling rate. For this test program the sampling rate is 65,536 samples/s. Consequently, 
the value is known within 0.015 ms (ss 1000/65, 536 ms). 

In order to estimate either time delay, D k g, a pattern matching function based on a regression line 
fit to the measured cross-spectrum phase g will be used. However, the regression line fit is not forced 
to have a zero intercept because 1) the pressure sensors are not the same; 2) the pressure sensors are 
not phase matched; 3) the phase measured at / = 0 will not be used in the time delay calculations; and 
4) while a large number of independent samples are used and the MSC is greater than 0.01, the phase is 
calculated with uncertainty fluctuations. In addition, the MSC plots show the phase related to the indirect 
combustion noise is in the frequency range from 16-208Hz. 13 ' 28 Consequently, only phase function points in 
that frequency range are used. 

1. Estimation of slope of cross-spectra phase and cost function 

Instead of using a simple linear regression model with a zero intercept at / = 0, a simple two parameter 
linear regression model is used. The procedure will be described using the combustor microphone and the 
far-field microphone. 


fik m(f) — a + i^fc m f + e (19) 

The parameters are estimated using the least square method that minimizes the sum of squares of the vertical 
distances ( S = £; (^ m (/j) -a-0 km /j) 2 ) ). 38 ~ 41 Let Xj = f, and y s = 4>% m (fj) where the time delay is 

D k m = 0 k m/(2 7r) , a — (f o and e is an error due to uncertainty fluctuations. 

Then 
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(20) 
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The residual variance is 


( 21 ) 

(22) 

(23) 

(24) 

(25) 

(26) 


(27) 

(28) 


_2 _ 
m 


£ j (<t>k m(fj)~ a- Pkm fj) 2 ) 


nk m 2 


(29) 


The standard error of the measurement is simply cr*, m . 

The cost function used to adjust the time delay Dm to minimize the slope of the linear curve fit to 
the thirteen phase measurements of interest and reduce the standard error is 


Cq — (0k m T ^ &k m) (30) 

The parameter u is selected to be about I.OjE — 04 so that if a value of <7*, m is 20, the value of ak m will not 
have an impact on the cost function until the slope, 0k m , has been reduced to l.E — 03. In order to prevent 
the search procedure from diverging two additional cost functions were defined. 


C m ax — F 2 ( Dk rn F) max ^ if Dk m P F) max (31) 

Cmin = F 2 (^mtn — Dk m) if Dk m ^ D m i n (32) 

where D max is some convenient number like 100 ms and D m i n is some convenient number like 70 ms and F 2 
is has some value like 10 or 100. These values serve to reduce the search region as needed. 

The total cost function, C, is 


C — Cq + C m a X C n 


(33) 


2. Parameter evaluation procedure 

The search technique used in this study is described by Powell 42 and a Fortran computer code for this 
algorithm is given by Kuester. 43 The search program evaluates, a new cost function with each iterative 
selection of Dk m . The signal r*. is measured at the combustor entrance at the same time as the signal rp is 
measured at the turbine exit and as each microphone signal is measured The subroutine that calculates the 
cost function time shifts the signal from the far- field microphone, r^, by an amount nk m using Fortran code 
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Ln993: do j=0,Lx 

nti) = r k (j+n km ) 

enddo Ln993 

where if the dominant low frequency noise source is indirect combustion noise, we let n k m = Int ( r s x D k m ) 
for the combustor signal case and nt m = Int(r s x xDg m ) for the turbine signal case. The actual value of 
D( m depends on the acoustic propagation velocity and distance along the paths from the turbine to the 
nozzle and from the nozzle through the jet at the nozzle exit and then to a microphone on the ground. The 
actual value of D k m depends on D? m and the velocity of the gas flow in the combustor and the distance 
in the combustor from the flame zone to the turbine. The cross-spectrum between the signals r k and r m 
is calculated. Next, a new cost function is calculated and the iterative procedure continues until it has 
converged. 

The resulting phase plot of <p k m | zero s i D p e over the frequency range from 16Hz to 208 Hz where 
the MSC is > 0.1 after convergence will show a a set of points moving above and below a horizontal line in 
a random fashion due to uncertainties caused by noise where 


4>k m I zero slope ~ arc t an 


^(G SlS 3| zer0 sloped) ex P(~^s lS 3 I zero slope (/))) 
s j 0 p e (/) exp (— i <Asis3| Z ero slope (-0)) 


(34) 


IV. Results 

A. Post-combustion residence time in the combustor 

The post-combustion (post-flame) residence times in the combustor, D k as determined for the 48 and 54 
per cent of maximum power operating condition cases is shown in Tables 2 and 3. These tables shows the 
propagation time for the turbine exit sensors is in the range of 86 to 94 ms and the indirect combustion noise 
signal takes 3.6 to 4.03 ms longer. Each table shows the actual number of counts and time in ms required 
to shift the far-held pressure signal to make it align with the combustion sensor signal. 

B. Analysis of distribution of post-combustion residence times in the combustor 

The distribution of the estimated post-combustion residence times in the combustor were analyzed using 
histogram density plots using the methods described by Fox 44 and by Fox and Weisberg. 45 The histogram 
is constructed so that the areas of the histogram bars sum to unity. Kernel density estimates are shown for 
the normal bandwidth (the line is of double thickness) and half the normal bandwidth (the line is of normal 
thickness). Also shown at the bottom of the plot are the location of the actual estimates as a scatterplot (or 
” rugplot” ) using a vertical bar | as the plotting symbol. A histogram of the post combustion residence times 
in the combustor at the 48 percent maximum power case is shown in Fig 2. Also shown is the mean value 
of the post-combustion time delay in the combustor at r mean = 3.960 ms. Using an alternative procedure 
with signals from only the combustor probe and turbine exit probes Miles 46 calculated a time delay of 
' r ciPi,T55i — 3.937 ms and a time delay of Tc 7 Pi,T 552 = 4.028 ms. These are also shown as vertical lines and 
are in close agreement with the values calculated using far-held microphones. 

A histogram of the post combustion residence times in the combustor at the 54 percent maximum 
power case is shown in Fig 3. Also shown is the mean value of the post-combustion time delay in the 
combustor at r mean = 3.791 ms. Using the alternative procedure with signals from only the combustor 
probe and the turbine exit probes Miles 46 calculated a time delay of rc/pi,T 55 i = 3.753 ms and a time 
delay of tc/pi,T 552 = 3.921 ms. These are also shown as vertical lines and are in close agreement with the 
values calculated using far- held microphones. The values obtained by Miles 13 15 ’ 28 ’ 47 using other methods 
are comparable to these. 

C. Phase and MSC measurements at 48 percent maximum power 

Aligned MSC and cross-spectrum phase measurements between the combustor pressure sensor (CIP1) and 
various far- held microphones are shown in hgures 4 and 5. Aligned MSC and cross-spectra phase measure- 
ments between turbine exit pressure sensor (T551) and various far-held microphones are shown in hgures 
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6 and 7. Aligned MSC and cross-spectra phase measurements between turbine exit pressure sensor (T552) 
and various far- field microphones are shown in figures 8 and 9. 

The magnitude-squared coherence function (MSC) plots will be discussed using logarithmic scales. 
Using logarithmic scales makes it easier to show that the analysis is based on MSC values above the threshold 
for unaligned signals. In addition to the MSC, also shown in Figs. 5, 7, 9, and 10 are MSC thresholds 
calculated from Eq (18). 

The 95 percent threshold confidence interval based on n s = 1120 independent samples is 2.67 x 10~ 3 . 
However, the spectra are calculated using a 50 percent overlap to reduce the variance and the 95 percent 
threshold confidence interval based on n s — 2240 samples is 1.337 x 10 -3 . These thresholds are shown in 
Figs. 4, 7, 9, and 10. These indicators show the MSC function is reliable up to about 400 Hz even though the 
pattern matching method used only cross-spectrum phase information from 16 to 208 Hz. However, MSC 
function is above 0.1 in a region from 16 — 208 Hz. 

The slope change in Fig. 4 above 200 Hz indicates that indirect combustion noise is important below 200 
Hz and direct combustion noise is important above 200 Hz. 

D. Core noise source 

In order to evaluate the relative strength of indirect combustion noise and the turbine exit coherent broadband 
noise which contains both indirect and direct combustion noise a comparison of aligned MSC between 
combustor pressure sensor (CIP1) and far- field microphone at 130 degrees with aligned MSC between turbine 
pressure sensors (T551) and (T552) and far-field microphones at 48 percent maximum power is shown in 
Fig 10 using both linear and logarithmic scales. Fig 10 b also shows MSC threshold plots for n s = 1120 and 
n s = 2240 samples. Figure 10 shows both indirect and direct combustion noise are important below 400 Hz. 

V. Concluding remarks 

A phase-based adaptive estimation pattern matching approach has been developed to estimate the 
post combustion (post-flame) time delay in the combustor from cross-spectral phase measurements made us- 
ing a combustor pressure sensor, two turbine exit pressure sensors , and far-field microphones by estimating 
the propagation time to the far field for the indirect combustion noise and the turbine exit coherent broad- 
band noise and taking the difference. The procedure accurately estimates the post-combustion/post-fiame 
residence time in the combustor. It is the presence of indirect combustion noise that makes the procedure 
feasible. While the method has general application, the method discussed works well in this particular case 
due to the nature of the data set. The number of independent averages is chosen to be large (n 8 ^1120) 
and the magnitude-squared coherence (MSC) is over 0.1 up to 208 Hz making the standard deviation of 
the cross-spectra phase small below 200 Hz. An accurate algorithm for calculating the post-combustion 
(post-flame) residence time in the combustor might be important in understanding the formation of thermal 
NO x in the combustor and verifying new and innovative combustor designs based on computational tools. 
Calculating the aligned MSC is a necessary step in applying array processing methods to the MSC data set. 
In general, the adaptive estimation technique developed herein is an effective tool to determine time delays 
between broadband coherent signals in noisy environments such as exist in real aircraft engines. 

VI. Acknowledgments 

This work was carried out under the NASA Fundamental Aeronautics Program, Fixed Wing Project, 
Quiet Performance Subproject. 


10 of 19 


American Institute of Aeronautics and Astronautics 



References 


1 Williams, J. E. F. and Howe, M. S., “The generation of sound by density inhomogeneities in low Mach number 
nozzle flows, J. Fluid Mech ., Vol. 70 Part 3, 1970, pp. 605-622. 

2 Pickett, G. F., 1 Core Engine Noise Due to Temperature Fluctuating Through Turbine Blade Rows,” Tech. Rep. 
AIAA-75-528, A1AA, 1975. 

3 Marble.. F. E. and Candel, S. M., “Acoustic disturbance from gas non-uniformities convected through a nozzle, 
Journal of Sound and Vibration , Vol. 55 No. 2, 1977, pp. 225-243, 

4 Cumpsty, N. A. and Marble, F. E., The Interaction of Entropy Fluctuations With Turbine Blade Rows; A Mech- 
anism of Turbojet Noise.” Proc. R. Soc. Lond. A., Vol. 357, 1977, pp. 323-344. 

5 Cumpsty, N. A. and Marble, F,, “Core Noise from Gas Turbine Exhausts,' Journal of Sound and Vibration , Vol. 54 
No. 2, 1977, pp. 297-309. 

6 Cumpsty, N. A., 'Jet Engine Combustion Noise: Pressure, Entropy and Vorticity Perturbations Produced by 
Unsteady Combustion or Heat Addition, Journal of Sound and Vibration, Vol. 66 No. 4, 1979, pp. 527-544. 

'Gliebe, P., Mani, R., Shin, H., Mitchell, B., Ashford, G., Salamah, S., and Connell, S., “Acoustic Prediction Codes," 
Tech. Rep. NASA CR-2000-210244, R99AEB169, General Electric Aircraft Engines, August 2000. 

s Mani, R., “Issues in Combustion Noise,” Annual Research Briefs 2007, Center for Turbulence 
Research, Department of Mechanical Engineering , Stanford University Stanford, Ca, 94305, US A, 2007, 

http: //ctr . Stanford. edu/ResBriefs07/22_mani_pp255_264 .pdf 

9 Bodony, D. J., “Scattering of an entropy disturbance into sound by a symmetric thin body," Physics of Fluids, 
Vol. 21 , 2009, pp. 096101-1-10. 

10 Leyko, M., Nicoud, F., and Poinsot, T., “Comparison of Direct and Indirect Combustion Noise Mechanisms in a 
Model Combustor,’’ AIAA Journal, Vol. No., November 2009, pp. 2709-2716. 

n Miles, J. H., VVasserbauer, C. A., and Krejsa, E. A., “Cross Spectra Between Temperature and Pressure in a 
Constant Area Duct Downstream of a Combustor,” AIAA-83-0762, NASA TM-83351, 1983. 

12 Miles, J. H., Wasserbauer, C. A., and Krejsa, E., “Cross spectra between pressure and temperature in a constant- 
area duct downstream of a hydrogen-fueled combustor,” Tech. Rep. NASA-TM-83463, NASA Glenn Research Center, 1983. 

13 Miles, J. H., Time Delay Analysis of Turbofan Engine Direct and Indirect Combustion Noise Sources, Journal of 
Propulsion and Power, Vol. 25 No. 1, January-February 2009, pp. 218-227. 

14 Miles, J. H., Core Noise Diagnostics of Turbofan Engine Noise Using Correlation and Coherence Functions,” 
Journal of Propulsion and Power, Vol. 26 No. 2, March/April 2010, pp. 303-316. 

15 Miles, J. H., “Separating Direct and Indirect Turbofan Engine Combustion Noise Using the Correlation Function," 
Journal of Propulsion and Power , Vol. 26 No. 5, September— October 2010, pp. 1144-1152. 

16 Weir(editor), D. S., “Engine Validation of Noise & Emission Reduction Technology Phase 1,” Tech. Rep. NASA/CR- 
2008-215225, Honeywell, May 2008. 

1 1 Etter, D. M. and Stearns, S. D., “Adaptive Estimation of Time Delays in Sampled Data Systems,” IEEE Transac- 
tions on Acoustics, Speech, and Signal Processing, Vol. ASSP-29 3, June 1981, pp. 582-587. 

18 Carter, G. C., “Time Delay Estimation for Passive Sonar Signal Processing,” IEEE Transactions on Acoustics, 
Speech, and Signal Processing , Vol. ASSP-29 3, June 1981, pp. 463-470. 

19 Knapp, C. H. and Carter, G. C., The generalized correlaton method for estimation of time delay,' 1 IEEE Trans. 
Acoust, Speech, Signal Processing, Vol. ASSP-24 No. 4, 1976, pp. 320-327. 

20 Azaria, M. and Hertz, D., “Time Delay Estimation by Generalized Cross Correlation Methods,” IEEE Transactions 
on Acoustics, Speech, and Signal Processing, Vol. 32 2, April 1984, pp. 280-285. 

21 Jacovitti, G. and Scarano, G., “Discrete Time Techniques for Time Dealy Estimation, -1 IEEE Transactions on 
signal processing , Vol. 41 No. 2, February 1993, pp. 525—533. 

22 Viola, F. and Walker, W. F., “A Comparison of the Performance of Time-Delay Estimators in Medical Ultrasound," 
IEEE '. Transactions on Ultrasonics, Ferroelctrics, and Frequency Control, Vol. 50 4, April 2003, pp. 392-401. 

23 Viola, F. and Walker, W. F., “A Spline-Based Algorithm for Continuous Time- Delay Estimation Using Sampled 
Data,” IEEE Transactions on Ultrasonics, Ferroelctrics, and Frequency Control, Vol. 52 1, January 2005, pp. 80-93. 

24 Piersol, A., “Time Delay Estimation Using Phase Data,' IEEE Transactions on Acoustics, Speech, and Signal 
Processing, Vol. ASSP-29 No. 3, June 1981, pp. 471-477. 

25 Chan, Y. T., Hattin, R., and J.B. Plant, ‘The least squares estimation of time delay and its use in signal detection,” 
IEEE Trans. Acoust. Speech, Signal Processing, Vol. ASSP-26 No. 3, June 1978, pp. 217-222. 

26 Zhen, Z. and Zi-qiang, H., The Geneaalized Phase Spectrtm Method for Time Delay Estimation,” Prloc. 
ICASSP’84,vol. 3, pp46. 2/1-4, March, 1984. 

27 Daneti, M., “On using phase data information for pipeline leak location,” Electrical and Electronics Engineers in 
Israel, 2008. IEEEI 2008. IEEE 25th Convention of, IEEE, 2008, pp. 494-498. 

28 Miles, J. H., “Spectral Separation of the Turbofan Engine Coherent Combustion Noise Component,” AIAA-2008-50, 
January 2008, NASA TM— 2008-215157-0,2008. 

29 Welch, P. D., “The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time 
Averaging Over Short, Modified Periodograms," IEEE Transactions on audio and electroacoustics, Vol. AU-15 No. 2, June 
1967, pp. 70-73. 

30 Stearns, S. D. and David, R. A., Signal Processing Algorithms Using Fortran and C, Prentice-Hall, Inc., 1993. 

31 Porat, B., Digital Processing of Random Signals, Dover Publications, NY, 1994,2008. 

32 Hayes, M. H., Statistical Digital Signal Processing and Modeling, John Wiley & Sons, 1996. 


11 of 19 


American Institute of Aeronautics and Astronautics 



33 Manolakis, D. G., Ingle, V. K., and Kogon, S. M., Statistical and Adaptive Signal Processing , Artech House, Boston, 
Massachusetts and London, United Kingdom, 2005. 

34 Bendat, J. S., “Statistical errors in measurement of coherence functions and input/output quantities,” Journal of 
Sound and Vibration , Vol. 59 No. 3, 1978, pp. 405-421. 

35 Bendat, J. S. and Piersol, A. G., Engineering Applications of Correlation and Spectral Analysis, John Wiley & 
Sons, 1980, pages 175, 167-172,215. 

36 Miles, J. H., “Aligned and Unaligned Coherence: A New Diagnostic Tool,” Tech. Rep. AIAA-2006-0010, AIAA, 
2006, Presented at the 44th AIAA Aerospace Science Meeting, 9-12 Jan 2006 Reno Hilton Reno, Nevada, also NASA/TM- 
2006-214112. 

37 Miles, J. H., “Estimation of signal coherence threshold and concealed spectral lines applied to detection of turbofan 
engine combustion noise, The Journal of the Acoustics Society of America, Vol. 129 No. 5, May 2011, pp. 3068-3081. 

38 Gibra, I. N., Probability and Statistical Inference for Scientists and Engineers, Prencice-Hall, Inc., 1973, page 

412-415. 

• 9 Edwards, A. L., An Introduction to Linear Regression and Correlation, W. H. Freeman and Company, 1976, pages 

24-31. 

40 Wackerly, D. D. and III, W. M., Mathematical Statistics with Applicatioins , Duxbury Pacific Grove, CA, 2002, 
pages 537-540. 

41 Chatterjee, S. and Hadi, A. S., Regressiion Analysis by Example (Fourth Edition), John Wiley & Sons, 2006, pages 

29-30. 

42 Powell, M. J. D., “An Efficient Method of Finding the Minimum of a Function of Several Variables without 
Calculating Derivatives, Computer Journal, Vol. 7 No. 2, July 1964, pp. 155-162. 

43 Kuester, J. L. and Mize, J. H., Optimization Techniques with Fortran, McGraw-Hill, 1973. 

44 Fox, J., An R and S-Plus Companion to Applied Regression, Sage Publications, 2002, page 86-89. 

45 Fox, J. and Weisberg, S., An R Companion to Applied Regression, Sage Publications, 2011, page 107-111. 

46 Miles, J. H., unpublished. 

47 Miles, J. H., “Separating Direct and Indirect Turbofan Engine Combustion Noise While Estimating the Post- 
Combustion (Post-Flame) Residence Time Using the Correlation Function,” Tech. Rep. NASA-TM-201 1-216248, NASA, 2011. 


12 of 19 


American Institute of Aeronautics and Astronautics 



Table 1. Spectral estimate parameters. 


Parameter Value 

Segment length, (data points per segment), N 4 096 

Sample rate, r s , samples/s 65 536 

Segment length, Tj = N/r s , s 0.062500 

Sampling interval. At = l/r s , s 1/65 536 

Bandwidth resolution, B e = A / = l/T^ = r s /N, Hz 16.0 

Upper frequency limit, f u = l/2Af = r s / 2, Hz 32 768 

Number of frequencies, Ly — / C /A/ = N/2 2 048 

Number of independent samples, n 3 1120 

Overlap 0.50 

Sample length, T tota i, s « 70 


Table 2. Post-combustion (post-flame) residence time in the combustor calculated from difference of propaga- 
tion time to far-field microphones of indirect combustor noise signal and acoustic signal at 48 percent maximum 
power. 


Far-held 
Microphone 
at 100 ft 

CIP1 

counts 

T551 

counts 

T552 

counts 

aGIPI 

^T551 

counts 

aCIPI 

552 

counts 

n T551 

U k e 

x 1000/2 16 
ms 

n T5 52 
U k t 

xlOOO/2 16 

ms 

O 

o 

Cfc 

6349 

6099 

6091 

250 

258 

3.814697 

3.936768 

110° 

6098 

5842 

5836 

256 

262 

3.906250 

3.997803 

130° 

5878 

5619 

5611 

259 

267 

3.952026 

4.074097 

160° 

5609 

5349 

5345 

260 

264 

3.967285 

4.028320 

mean (r) 

3.960 ms 







standard deviation (t) 

0.079 ms 








13 of 19 


American Institute of Aeronautics and Astronautics 





Table 3. Post-combustion (post-flame) residence time in the combustor calculated from difference of propaga- 
tion time to far-field microphones of indirect combustor noise signal and acoustic signal at 54 percent maximum 
power. 


Far-field 
Microphone 
at 100 ft 

CIP1 

counts 

T551 

counts 

T552 

counts 

^T551 

counts 

\C1P1 

552 

counts 

r>T551 
n k l 

xlOOO/2 16 

ms 

n T552 
U k l 

xlOOO/2 16 

ms 

90° 

6367 

6131 

6124 

236 

243 

3.601074 

3.707886 

O 

O 

6112 

5868 

5859 

244 

253 

3.723145 

3.860474 

130° 

5891 

5640 

5636 

251 

255 

3.829956 

3.890991 

160° 

5620 

5366 

5368 

254 

252 

3.875732 

3.845215 

mean (r) 

3.792 ms 







standard deviation (r) 

0.103 ms 








Sensors in aft fan 
bypass access panels 


High-temperature sensors 
with air cooling at turbine exit 



-temperature sensor 
with air cooling in 
combustor igniter port 


Figure 1. Honeywell TECH977 engine. 
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time delay, ms 


Figure 2. Histogram of post-combustion residence times in combustor in ms at 48 percent maximum power. 
Shown at the bottom of the plot are the location of the actual estimates as a scatter plot using a vertical 
bar | as the plotting symbol. Also shown are r me0 n (blue dashed), t C ipi,t 551 (red dashed), and t C ipi,t 552 (red 
dashed) . 



Figure 3. Histogram of post-combustion residence times in combustor in ms at 54 percent maximum power. 
Shown at the bottom of the plot are the location of the actual estimates as a scatter plot using a vertical 
bar | as the plotting symbol. Also shown are T mean (blue dashed), t C ipi,t551 (red dashed), and t C ipi,T5S2 (red 
dashed). 
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Frequency, Hz 


Figure 4. Aligned cross-spectra phase angle between combustor pressure sensor (CIP1) and far-field micro- 
phones at 48 percent maximum power. 



Figure 5. Aligned coherence between combustor pressure sensor (CIP1) and far-field microphones at 48 percent 
maximum power. 



Figure 6. Aligned cross-spectra phase angle between turbine pressure sensor (T551) and far-field microphones 
at 48 percent maximum power. 
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Figure 7. Aligned coherence between turbine pressure sensor (T551) and far-held microphones at 48 percent 
maximum power. 



Figure 8. Aligned cross-spectra phase angle between turbine pressure sensor (T552) and far-field microphones 
at 48 percent maximum power. 



Figure 9. Aligned coherence between combustor pressure sensor (T552) and far-held microphones at 48 percent 
maximum power. 
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(a) Linear scale. 



(b) Logarithmic scale 


Figure 10. Comparison of aligned coherence between combustor pressure sensor (9) and far-field microphone 
at 130 degrees with aligned coherence between turbine pressure sensors (10) and (11) and far-field microphones 
at 48 percent maximum power. 
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